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He-Ne ring laser gyroscopes are, at present, the most precise devices for 
absolute angular velocity measurements. Limitations to their performance 
come from the non-linear dynamics of the laser. Following the Lamb semi- 
classical theory, we find a set of critical parameters affecting the time stability 
of the system. We propose a method for estimating the long term drift of 
the laser parameters and for filtering out the laser dynamics effects from the 
rotation measurement. The parameter estimation procedure, based on the 
perturbative solutions of the laser dynamics, allow us to apply Kalman Filter 
theory for the estimation of the angular velocity. Results of a comprehensive 
Monte Carlo simulation and results of a preliminary analysis on experimental 
data from the ring laser prototype G-Pisa are shown and discussed, (c) 2012 
Optical Society of America 



1. Introduction 

Ring laser "gyros" are standard sensors in estimating rotation rates relative to the inertial 
frame, with many applications ranging from inertial guidance [l], to angle metrology p], 
geodesy [3111], geophysics [5], as well to fundamental physics [5l|6]. In the near future their 
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application is foreseen to improve the performances of advanced gravitational waves detectors 
and also to provide ground based tests of General Relativity [8]. 

In a ring laser two oppositely traveling optical waves resonate inside a polygonal closed 
path. In an inertial frame, each beam follows a path of the same length. When the cavity 
rotates, a round-trip time difference between the co-rotating and counter-rotating beams 
occurs, as they experience a longer and a shorter path, respectively. This translates into a 
frequency difference of the two beams. Such frequency difference contains the information 
about the rotation rate of the reference frame. 

This physical phenomenon is known as Sagnac effect and the Sagnac frequency (i.e. the 
frequency of the beat signal between the two beams), reads 



where A and L are the area and the perimeter of the cavity respectively, A is the wavelength 
of the laser beam, n is the normal vector of the plane of the ring cavity, Q, is the rotation 
rate vector. In a real ring laser the laser dynamics are source of systematic errors in the 
estimate of Ug. The laser dynamics are determined by a set of non linear equations that de- 
pend on parameters that are slowly varying, following random changes of the environmental 
conditions (mainly temperature and atmospheric pressure). In addition, the control of some 
parameters (e.g. laser emission frequency or light intensity) in a closed loop way may en- 
hance the drifts of the others. Thus the simple Eq. ([1]) must be modified to take into account 
such complex behaviors (see for example |5]) which result in i) corrections to scale factor 
due to fluctuations of the laser gain, cavity losses and frequency detuning; ii) null shift due 
to any cavity non- reciprocity; and in) non linear coupling between the two laser beams due 
to backscattering (at present, the most important instability source). Therefore, the beating 
signal can be represented as the sum of Ug , white noise and terms arising from the non linear 
dynamics. Usually, the last ones are removed by long term correlations with a set of auxil- 
iary sensors information, and by modeling the ring laser behavior with respect to the ideal 
case [9]. However, the non linearity of the dynamics limits the effectiveness of this approach. 
In this paper we propose to increase the ultimate resolution and the time stability of ring 
lasers by the estimation of the parameters of their system, and the subsequent application 
of an Extended Kalman filter (EKF) [lO]. 

The proposed method has a very wide range of application. In fact, the parameter estima- 
tion and the dynamical filtering could improve the response of large size ring lasers, which 
have already reached accuracy of ~ 1 part in 10^, and so providing unique informations for 
geophysics and geodesy. Middle size rings, with sides of ~ 1 m, which are more affected 
by backscattering, will be improved as well. Such instruments are more suitable for geo- 
physics apphcations (i.e. rotational seismology), and for application to gravitational waves 
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interferometers (e.g. local tilts measurements). 

The paper is organized as follows. In Section [2] we discuss the non-linear dynamics of 
a He-Ne inhomogenously broadened ring laser. Here the critical parameters affecting ring 
behavior are presented. Section [3.Bl describes the implementation of the parameter estimation 
procedure and the Extended Kalman Filter (EKF) algorithm [10]. In Section H] we presents 
the results of the parameter estimation procedure, and the performance of EKF for rotation 
estimation. The procedure has been tested on the experimental data of G-Pisa, and results 
are reported in Section [5l Finally, our conclusions are drawn in Section [6l 

2. Dynamics of a Ring Laser 

The differential equations of the ring laser dynamics were first derived long ago, using a low 
gain self consistent model by E. Lamb [11] for linear laser, and then extended to rings by F. 
Aronowitz [12] 
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where Ii^2, ^iid tp are the dimensionless light intensities, the instantaneous phase difference, 
and the instantaneous circular beat frequency of the counter-propagating waves, respectively. 
Here Ug = 2'kVs is the rotation rate in Eq.([T]), and 011^2-, <^i,2i Pi.2^di2,2ii'Ti2,2ii^i,2i ^ are the 
Lamb parameters. It is worth mentioning that ai^2, (^1.2 are the amplification minus losses, /3i.2 
is the self saturation, ^12.215 ''"12,21 describe cross-(mutual-)saturation, ri 2, s are the amplitude 
and the relative phase of the backscattered waves, respectively. A more detailed explanation 
of the physical meaning of these parameters can be found in Appendix [X] 

2. A. Lamb parameters effects on the gyroscope performances 

The study of Eqs.([2]) has been conducted by several authors in the past [T314T7] . both with 
numerical and analytical approaches. The analytical solution cannot be found in the most 
general case, but only under certain approximations about the reciprocity of the system. 
Approximated analytical expressions for the time evolution of the Sagnac phase provide, 
nevertheless, an useful reference to better understand the role of the Lamb parameters noise 
on the estimation of the angular velocity Us- 

We present in the following the periodic solution of Eqs.([5]) in the case where: /1//2 = k 



3 



and Ti2 = T21 = 0. In this case the phase equation takes the form: 



ip = Us - Y 



Tik sin(?/; — e) + — sin(?/; + e) 



(3) 



and admits the following solution: 



ip{t) = 2 arctan 



VLli + rip tan {\^pt) 



(4) 



where VL^i = c/L{kri + r2/A;) cos(£), = c/L{r2/k — /cri) sin(e) and Qp = 



where oo denotes the detected Sagnac frequency and oobs = —^L2COSUJst — Q^LiStinoOgt 
represents the frequency modulation of the Sagnac signal. 

In Figured] we report the results of a Monte Carlo simulation of lO^s of ring laser dynamics 
evolution with /1//2 = k. We considered the following noise sources in the system: a white 
frequency noise with standard deviation of 10^^ on u)(t) (mimic of the output of AR(2) 
frequency detection algorithm [20]) and a random- walk noise on the parameters ri 2, k, and 
e. The Allan variance of uj{t) has been calculated for four different cases, denoted with (a), 
(b), (c) and (d). 

In case (a) (ri 2,A;,£:) are independent and vary in random walk with a step size of (2 ■ 
10~^, 10"'^, 10~^) respectively. In case (b) the only varying parameter is k, with a step size 
of 10~^. In case (c) all parameters vary as in (a), but the processes ri and r2 have been 
correlated with a correlation coefficient of 0.9 while e varies around the nominal value of 
Orad with a random walk step size of 10~"^rad. In case (d) all parameters vary as in (b), but 
the process e varies around the nominal value of 7r/2rad with a random walk step size of 
10"^ rad. 

It can be easily observed that the noise contribution coming from the parameters fluctua- 
tion is transferred to the noise of the measured Sagnac frequency exhibiting the same random 
walk plus white noise pattern. The relative noise on the laser parameters is converted into 
frequency noise by the factor c/L meaning that the larger is the cavity perimeter, the larger 
is the rejection of the laser parameters noise. In addition, it is worth noticing that the 
backscattering phase e plays a crucial role in transfering the fluctuations of ri,r2 on u. It 
determines a strong reduction of the output noise for values close to e = tx/2 (trace (d)). 
In this regime, also known as 'conservative coupling regime' [H], the backscattered pho- 
tons interact destructively and their influence on the nonlinear interaction between the two 
intracavity beams and the active medium is minimized. 




uj{t) 



~ cjg — VtL2 COS ujgt — sinwgt 



(5) 
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10^ 10^ 10^ 10* 

Fig. 1. Simulated Allan deviations of the estimated rotation rate. See the text 
for details. 

In the next sections we will show that a statistical filtering procedure is able to identify 
the Lamb parameters, and to remove their slow drifts by using extended Kalman filtering. 
Thus the maximum resolution (i.e. the minimum value of cr^(T)) and the time stability (i.e. 
the value of r where the minimum is attained) of a ring laser can be significantly improved. 
In fact, the instantaneous Sagnac frequency will depend in general on the full system state 
[/i(t), l2{t), ip{t)], and so dynamic Kalman filtering can be more effective in estimating Ug 
than other approaches that rely on ilj{t) only (e. g. the standard AR(2) method for frequency 
estimation |20|). 

3. Dynamics of G-Pisa 

G-Pisa is a prototype middle size He-Ne ring-laser. The main characteristics of its optical 
cavity are reported in table [1] while its experimental setup is sketched in figure [2l 

The laser operation is controlled by two feedback loop systems, dedicated to the active 
stabilization of the optical frequency and the optical power [6]. The first loop keeps constant 
the frequency difference between the ring laser clockwise beam and a reference laser. It 
acts on two piezoelectric transducers (PZTi 2) moving two opposite cavity mirrors along the 
cavity diagonal. The second loop regulates the RF discharge power in order to keep constant 
the power of the clockwise ring laser output. 

G-Pisa dynamics can be derived from Eqs.(l2]) taking into account the spectroscopic prop- 
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Ii{n} 




IBS 



Fig. 2. G-Pisa experimental setup. The cavity vacuum chamber is entirely filled 
with a mixture of He-Ne and does not contain any intra-cavity element except 
for the four mirrors. Sju}: Sagnac interference signal; Ii{n}: counter-clockwise 
single beam intensity; l2{n}: clockwise single beam intensity; Ib: optical beat 
intensity; RFD: radio frequency discharge; IBS: intensity beam splitter; HWP: 
half wave-plate; PZT: piezoelectric transducer. 
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G-Pisa 


Geometry 


Cavity 


square 


Side length 


1.35 m 


Latitude 


43° 40' 35.86"N 


Cavity mirrors 


Radius of curvature 


4 m 


Total losses 


3.7 ppm 


Transmission 


0.25 ppm 


Scatter+absorption 


3.5 ppm 


Optical properties 


Wavelength 


632.8 nm 


Output power 


1.6 nW (single mode) 


Spatial mode 


TEMoo 


Beam waist (s,h) 


(1.97 mm, 2.43 mm) 



Table 1. Main nominal characteristics of the "G-Pisa" apparatus. Mirrors char- 
acteristics are the nominal ones (manufacturer information) and refer to the 
single mirror for an incidence angle of 45°, and s-polarized light at 632.8 nm. 
Beam waists are reported as 4 times the 1/e^ beam intensity radius. 



erties of its active medium and the constraints imposed by the stabilization loops. The use of 
a special gas mixture, containing ^''A^e and ^^A^e isotopes at 50:50 ratio, avoids mode com- 
petition effects, while the perimeter control, by keeping constant the laser optical frequency 
avoids mode-jumps. In the Appendix [A] we show that (3i = P2 = P, o"! = (^2, T12 = T21 = 0, 
and 6*12 = 6*21 = hold for the closed loop operation of G-Pisa. Thus the equations of the 
dynamics reduced to: 



L L 

c 

L L 



aJi - I3ll + 2r2 v7^cos(?^ + e) 
Oi2h - Pll + 2riA//i/2 cofi{ip - e) 



(6) 



ri 



sin(?/' — e) 



r2^/ — sin(?/' + £:) 



In Table [2] we report the typical Lamb parameters of G-Pisa, with a typical round-trip gain 
of G = 3 X 10~^. It is worth noticing that, as the perimeter is kept constant by moving the 
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ring mirrors, the backscattering phase angle e is not fixed, but it ranges from (conservative 
couphng) to 7r/2 (dissipative couphng) |14] . 



Parameter 


Value 


Us 


2n 107.3 Hz ~ 674.5 rad/s 


a 


4.3 X 10-9 


0.24 rad/s 




1.5 X 10-^ 


825 rad/s 


r 


2 X 10-^ 


11.1 rad/s 



Table 2. Typical values of Lamb parameters used in the simulations of G-Pisa 
dynamics with G ~ 3 ■ 10~^. In the first column the units are adimensional 
Lamb units. To compare the order of magnitude of effects with Ug, in the second 
column we report such parameters multiplied by the free spectral range c/L. 



3. A. Steady state approximated solutions 

To provide suitable algorithms for parameter estimation, we study the steady state regime 
of Eqs.(l6]). By inspection of the right hand side of Eqs. ([6]), one finds that the general steady 
state solutions are periodic. In particular, without backscattering (ri^2 = 0), Eqs. ([6]) exhibit 
steady state solutions of the type 



hit) = f 

m = f (7) 

1p(t) = Ust 

for t — J- OO. In the presence of backscattering, the above solutions switch to periodic steady 
state solutions and exhibit oscillatory behaviors, and the backscattering can be treated as 
a perturbative sinusoidal forcing term. We can study the system oscillation around its un- 
perturbed steady state by means of the time dependent perturbation theory To this 
aim we introduce the expansion parameter A, which is assumed to be of the same order of 
magnitude of ri_2 and write: 



Xi(A,t) 
'^(A,t) 



OO , ^ 



k=0 

OO 



A;! 



k=0 

OO 



k=0 



k\ 
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For k = 0, substituting the latter into the Eqs.([6]), we recover the solution ([7]) with the 
positions lf'\t) = l2'\t) = 02/ /3, ip^'^''{t) = Ugt. The approximated solutions can be 

calculated iteratively from the series expansion in powers of A of ([8]) into the dynamic of 
Eqs([6]). A second order approximations of the solutions reads: 



^1 W- ^ + 2r2v/«ia2 7 r 2 sm(2e) 

^ , , a2 , ^ . «2Cos(£-a;,t)-(^)sin(5-a;,t) rir2(c/L) . ,^ , 

2^2(t)^ -T + 2ri^aia2 7 ^ x + 2 sm(2e) 

/3 /3r«2^(-^12^ 



, 2nr2(c/L)2cos(2e), , f^cos{e - u,t) + r,^ f^cos{e + uj,t) 

(9) 

where we made the additional approximation of keeping the leading terms in w", with n <2. 
Solutions ([H]) show a correction to the mean intensity level and pushing and pulling in the 
phase difference, as well as the presence of the first harmonic of Us- 

3.B. Parameter Estimation and Kalman filtering 

The Sagnac phase can be conveniently estimated by means of the Hilbert Transform (HT) 
of the interferogram S{t) ~ sin?/'(t) which is routinely acquired during ring laser operation 
together with the monobeam intensities and l2{t) |22]- Since ring laser signals are 

sampled, we shift to the discrete time domain, and from here on we denote time dependent 
intensities and Sagnac phase as {J2(n)} and {n G N), respectively. 

Perturbative solutions and the least squares methods provide the main tools for estimation 
procedure of the parameters ai 2, ^"1.2 and e. For sake of clarity, we re-parametrize the 
measured intensities and Sagnac phase in the form 

Xi(t) = Ji + 2i sin (wt + 0i) 

J2(t) = /2 + ^2 sin (wt + 02) (10) 
^{t) = ut 

as the quantities I12, ii,2 and 0i 2 can be readily estimated from {Ii 2{n)} and {S'(n)}. 
In fact, the mean intensities Ji^2 can be estimated by computing the sample average of 
{Ji 2(n)}. Moreover, the averaged modulation amplitudes ii 2 and phase difference 0i — 02 
at the fundamental frequency can be estimated by means of a digital lock-in procedure, 
which calculates the "in-phase" and "in- quadrature" components of {/i 2(^)}. The reference 



9 



complex signal for the digital lock-in is given by the HT of {S{n)}. Averages are taken over 
a time interval where Lamb parameters remain fairly constant. The estimation procedure of 
Lamb parameters can be conveniently divided into two steps: 

1. The first step is to estimate the phases of the intensities. From Eqs.Q in the approx- 
imation j^YT >> «! 2 we have 

(c/L) 

Ui = e 
[02 = -e . 

Thus we can immediately identify the backscattering angle from 0i 2 as 



e = 




where the hat symbol ^ denotes an identified parameter. 

2. In the second step the remaining Lamb parameters are obtained by least squares meth- 
ods. In fact, starting from the periodic steady state solutions f llOp . we can form the 
squared residuals 



A(ai,a2,ri,r2) = — x 

-2^/1^ 



J r^~L r^^^ ~ ^^^^ ^ 2r2 VX1X2 cos{ut + ^)\ \ + 



X 



^2X2 - /3X2 + 2riA/XiX2 cos{ut - e) 



dt 



(12) 



averaged over a period 211 /uj. Minimization of A(ai^2, ''"2,1) yelds the best linear estimate 



of ai^2 and ri 2; from the conditions dA/dai 
dA/dri = 0, we get 



0, dA/da2 = 0, dA/dr2 = 0, and 



?2W 



2{c/L)^Ij2 
2(c/L)v^ 



ilt2UJ 

4(c/L)/2 

«1«2W 

4(c/L)Ji 



sin 2£ 
sin2e 



(13) 
(14) 
(15) 
(16) 
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which fulfill the parameter estimation procedure via the second order approximation. It is 
worth noticing that one can increase the precision of the identified parameters by evaluating 
solutions of Eqs. (2) of higher order in A and increase their accuracy by increasing the 
averaging time span. 

4. Simulation results for G-Pisa 

We briefiy describe the specific implementation of the parameter estimation procedure for the 

G-Pisa ring laser, where the data are acquired at a sampling frequency of 5 kHz (Tg = 200 fis). 

To remove the oscillating component, intensity signals are low-pass filtered with a first order 

Butterworth filter with 1 Hz cutoff frequency. The quantities I12 are estimated by averaging 

the decimated intensities over a time interval of 10 s (i.e. 5 x 10^ samples). On the other side, 

to calculate the modulation ii^2 and phases (f)i^2, the intensities are first band-passed around 

the fundamental Sagnac band [95 -r- 125] Hz by means of a Butterworth filter, and decimated 

by a factor 2. The decimation procedure has been carried out by the tail recursive routine 

"Zoom and Decimation of a factor 2"" {ZD{n)), where each iteration step is composed by a 

^3 _^ 22"^ + 2z + 2 

half band filter stage with discrete transfer function H(z) = , followed by 

a downsampling by 2. The ZD(n) procedure ensures a linear phase filter response at least 
for n = 3 iterations, as no appreciable phase distortion was observed in simulated sinusoidal 
signals. The resulting data are then demodulated with a digital lock-in using as reference 
signal the discrete HT of the interferogram, and setting the integration time to 10 s. A 
schematic of the parameter estimation procedure is reported in Fig. [3l In addition, the phase 
of the two monobeam oscillating components is determined by the discrete HT, and their 
difference is estimated by unwrapping the phase angle and taking its average over 10 s. As a 
concluding remark on the parameter estimation procedure, we mention that the problem of 
filtering very long time series, has been solved by the "overlap and save" method [21], which 
is an efficient algorithm for avoiding the boundary transients due to finite length of digital 
filters. 

Reliability of the parameter estimation routine is tested by Monte Carlo simulation of the 
dynamics of Eqs. ([6]) followed by the estimations of Lamb parameters from the simulated 
time series of Ji 2 and tp. We run lO'^ simulations of the dynamics of Eqs. ([6]) allowing ai 2 
and ri 2 to vary according to normal distributions with mean as in Tab. [2] and standard 
deviation equal to 10% of their means. In addition, /3 is assumed constant, and e uniformly 
distributed in [0, 7r/2). In each simulation, we have compared the numerical RK4 solution of 
Eqs. ([6]) and approximated analytical solution ([9]) evaluated with the same Lamb parameters. 
We found that they are in a very good agreement, with means of the relative errors on 
l2{t) and tlj{t) of —6.4 x 10"'', —6.2 x 10"'' and —1.5 x 10"^, and standard deviations 
of 4.6 X 10~^ 4.6 X 10~6 and 1.3 x 10"^ respectively. 
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Fig. 3. Schematic of the parameter estimation procedure, where LP — > low- 
pass Butterworth filter, BP — )■ bandpass Butterworth filter, ZD — )■ Zoom and 
Decimation routine, HT — )■ Hilbert transform (see text). 



To numerically assess the performance of the parameter estimation procedure, we run a 
simulation of 6 hours where ai 2, 2, e and Ug fluctuate following independent random walk 
processes with self-correlation time of 1 hour. To reproduce the experimental behavior of a 
ring laser, the time drift of Ug, which mimics the effects of local tilts and rotations, is a factor 
of 5 lower than the auto-correlation time of the other parameters. We superimposed to the 
simulated data an additive white noise, with SNR = 10^ for the beam intensities and SNR 
= 5 X 10'^ for the interferogram. Such order of magnitudes are routinely achieved in large 
ring laser [25] and in G-Pisa j6]. The results we got are summarized in Fig. HI Fig. [5] and 
Fig. El 

The overall accuracy of the Lamb parameter estimation procedure is good, with a relative 
standard deviation of 3 x 10~^ and 4 x 10~^ in the estimation of ai 2, and ri 2, respectively. 
The absolute error in the estimation of the backscattering phase is 3 x 10^'^ rad. The attained 
accuracy is not far from the lower bound ~ 10"^ associated to the level of the observation 
noise of {Ji(n)}, {J2(n)} and {ilj{n)}. 

4 -A. Estimation of Ug by Kalman filtering 

Knowledge of the Lamb parameters ai^2 ,^1,2 and e, together with the /3 parameter separately 
acquired, allow us to set up an Extended Kalman Filter [10] for the estimation of the rotation 
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100 



lOOr 




Fig. 4. Histograms of the relative errors (ai^2 — cn,2)/cn,2 that affect the esti- 
mation of Gain minus losses parameters calculated with 2 x 10^ realizations 
of the ring laser dynamics, (a) Histogram relative to cti : mean 1.4 x 10~^ and 
standard deviation 2.9 x 10"'^; (b) histogram relative to ai : mean —2.5 x 10~^ 
and standard deviation 3.9 x 10"'^. 




0.02 



Relative Error Relative Error 

(a) (b) 

Fig. 5. Histograms of the relative errors (fi 2 — ''"1,2) /'"i, 2 that affect the es- 
timation of backscattering coefficients calculated with 2 x 10^ realizations of 
the ring laser dynamics, (a) Histogram relative to r\ : mean 1.1 x 10"^ and 
standard deviation 4.6 x 10"'^; (b) histogram relative to r2 : mean 1.3 x 10"^ 
and standard deviation 3.2 x 10"^. 
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Fig. 6. Histogram of the absolute errors ^ — e that affect the estimation 
of backscattering phase calculated with 2 x 10^ realizations of the ring laser 
dynamics; mean —4.3 x 10~^ rad and standard deviation 2.8 x 10"'^ rad. 

rate Us- 

The EKF state variables are the W vector X(t) = [h{t), hit), ip{t)Y ■ The dynamics 
model is given by Eqs.dH]), with the addition of the model error as a zero mean, white, 
stochastic vector field v(t) with variance Var[v(t)] = Q, where Q is a 3 x 3 covariance matrix 
that accounts for the effects of unmodeled dynamics, for instance, identified parameter errors, 
calibration errors, and numerical integration inaccuracies. The EKF prediction step, which 
corresponds to the integration of Eqs.(l6]) over the time interval T,, is carried out using the 
RK4 Runge-Kutta routine. 

In the discrete time domain, the model of the measurement process reads {y{n)} = 
{X(n)} + {w(n)}, where w(n) is zero mean, white, stochastic vector field (observation 
noise) with variance \^ar[w(n)] = R, and i? is a 3 x 3 covariance matrix. In the standard 
experimental set up of ring lasers [Ii{t), hit), ipit)] are measured by independent sensors, 
and so we can assume that R is diagonal, with diagonal elements the observation noise 
variances cr|^, a]^, cr^ which can be conveniently calculated through the level of white noise 
in the power spectrum of {y(ri)}. 

The backscattering frequency is estimated from the filtered channels Ii^2{n), ip{n), the 
identified parameters 0-12, 2, ^, and the exogenous parameter /3 as 

^ c 

where, for simplicity, we have dropped the index (n) from time series. The Sagnac frequency 
is then estimated from the difference Qs = 4' ~ ^bs, where the numerical derivative of ip has 




fi W — sm{ip — e) +r2\ — sm{ip + e) 
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Fig. 7. Allan standard deviation of the rotation rate estimated by AR(2) 
method (circles) and Extended Kalman Filter (open diamonds) using 2 x 10^ 
seconds of the Monte Carlo simulation with random walk of Lamb Parameters. 
For comparison, we also plot the Allan standard deviation of the simulated ro- 
tational drift (open triangles), and the Allan standard deviation of the EKF 
estimation after the subtraction of the rotational drift (dash dotted line). 

been computed by the "5 point method" [H] designed to reject the derivative amplification 
of the noise. 

The capability of the EKF in increasing the time stability and the resolution of the gyro- 
scope has been tested with a 6 hours simulation of the ring laser dynamics with parameter 
variations as in the parameter estimation tests. The results are summarized in FigJT] where 
we compared the Allan variance of AR(2) and EKF frequency estimations. We conclude that, 
for this simulation with typical parameters of middle-size rings, the rotational resolution in- 
creases by a factor of 10 while the minimum of the Allan standard deviation shifts from 60 
s to 360 s. 

5. Application to real data 

To run the parameter estimation routine on the G-Pisa data, the photodetector signals V\;2,{t) 
have been converted to dimensionless intensities h,2{t) (Lamb units), using the following 
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relation 



amb~ 



^LambPout 1,2? 



phO-eff 

where: Gph = 10^ V/A is the photo-amphfiers gain, Oe// = 0.4A/W is the quantum efficiency 
of the photo-diodes, and CLamb = 3.5 x 10^ is the cahbration constant to Lamb units ( 
see Appendix Ej); Pout 1,2 is the output power in Watt. The parameter estimation for G-Pisa 
is completed by the acquisition and calibration of the /3 parameter. 

In Fig. |8] we show the comparison between the real time series measured on the G-Pisa 
ring laser and the calculated signal by the model after the parameter estimation according 
to the scheme in Fig. [3l 



Monobeam Intensity 1^ 




0,03 0,04 
Monobeam Intensity I 




0,02 



0,03 0,04 
Sagnac Interferogram 



0,06 




Fig. 8. Extended Kalman Filter estimation of l2{'n) and sinipln). Circles: 

G-Pisa raw data sampled at 5 kHz. Continuous line: filter output. 



In Fig. |9]we show the time series of the identified parameters for the G-Pisa ring laser and 
the calibrated parameter /3, using the routine described in Section 13. B[ 

After the estimation of the G-Pisa parameters, we apply the EKF to the light intensities 
and to the interferogram (see |22]). However, the implementation of the EKF requires an 
estimation of the covariance matrices Q and R of observation and model errors. Typically Q 
and R are considered as tuning parameters and set on the base of trial-and-error procedures. 
In fact, we started from an initial raw estimation for the diagonal elements of Q and R using 
simulations and power spectra of Ji 2 and ip, respectively. Then we tuned these values search- 
ing for the minimum of the Allan variance of Ug and came to Q = diag(10~^, 10"^, 10"^'') and 
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Fig. 9. Time series of ai^2, fi,2 and e and of /3, respectively estimated and 
calibrated using 6 hours of experimental data of G-Pisa. 



R = diag(10~*, 10~*, 10~^). The performance of parameter estimation and EKF were limited 
by the environmental conditions of G-Pisa, e.g. local tilts and spurious rotations induced by 
the granite slab that support the instrument and some electronic disturbances, as it can be 
seen in Fig. [TOl where we report the power spectrum of {S{n)}. 

In figure [11] we report the Allan standard deviation of the Sagnac frequency estimated 
with AR(2) and EKF. An increase of a factor of 1.5 in rotation-rate resolution and of a 
factor of 2 in the time stability is observed. 

6. Conclusions 

A full model of the ring laser dynamics has been studied and applied to the estimation and 
removal of the long term drift in the laser parameters. The proposed data processing tech- 
nique is based on the estimation of the parameters appearing in the Lamb equations for a 
He-(^°Ne Ne) ring laser. Results of Monte Carlo simulations supported the viability of 
the parameter estimation, and yielded a relative estimation error of the order of 3 x 10"'^ 
for ai^2 and ri^2 and an estimation error ~ 10^^ rad for e. The accurate estimation of the 
Lamb parameters allows for the application of the Kalman Filter for the estimation of the 
rotation rate from the Sagnac frequency. Simulations showed a significant improvement in 
the frequency estimation by EKF compared to AR(2) method. Preliminary results on data 
from the ring laser prototype G-Pisa, presently strongly affected by local rotational noise. 



17 




80 90 100 110 120 130 140 

Frequency [Hz] 



Fig. 10. Power spectrum of the interferogram data around the Sagnac fre- 
quency ~ 107.3 Hz. 

make us confident about the rehability of our approach in the presence of unmodeled exper- 
imental noise and cahbration errors. Our approach can be further improved, and different 
mathematical tools can be used. For instance, the numerical integration method RK-4 can be 
substituted by geometrical integrators, or modified in conservative routines [27]. EKF can be 
also improved by increasing the state dimension, modifying the observation model and the 
estimation of Ug. As a final remark, we note that ring lasers achieved world record in rotation 
sensitivity, accuracy and time stability with sophisticated hardware, accurate selection of the 
working point of the He-Ne laser, despite a very basic off-line analysis. In this paper, we have 
shown that the parameters of the ring laser dynamics can be identified and their effects on 
resolution and time stability removed notwithstanding the system non linearities. We think 
that data analysis will cooperate more and more with ring laser hardware in pushing the 
resolution and the time stability of ring lasers beyond the current limitations. 
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A. Lamb Coefficients calculation 



The two intensities Ji and I2 in Eqs.([2]) are expressed in Lamb Units 

\^^ab\'^{la + lb) Pout 1,2 



J _ l/iq^l {la + lb) ^2 
^ri^lalblab 



4:h'^lalblab 2ceoSbT 



CLarnbPout 1,2 



where 



is the electric-dipole matrix element between the laser states a = and b = 2p^ (i.e. the 
upper and the lower of the laser energy levels), 7^ and 7b are the decay rates in Paschen 
notation, 77 = '~fab/'^D is the ratio between the homogeneous and the Doppler broadening of 
the laser transition. Here Yd = ^y2khTp/mne/^, where ki, is the Boltzmann constant, m„e 
is the atomic mass of Neon and Tp is the plasma temperature. Aik is the radiative decay 
rate between the laser levels, Ei^2 are the electric fields amplitudes. Pout 1,2 are the output 
powers, eo is the dielectric constant of vacuum, Sb is the area of the transverse section of the 
beam, T is the transmission coefficient of the mirror, and h is the reduced Plank constant. 
The table [3] contains the reference values of the above quantities for a Doppler broadened 
active medium in presence of collision, according to refs. [11] and |26] . 



Sb 


18 ■ 10-6 m^ 


p 

out 


1 - 10 nW 


P 


5.25 mbar 


la 


12 MHz 


lb 


127 MHz 


lab 


234 MHz 


T 
j-p 


450 K 


IJ-ab 


3.2 ■ 10-30Cm 



Table 3. The G-Pisa laser parameters. 



The coefficients in Eq.([2]) can be calculated by means of the plasma dispersion function, 
which is the function associated to the broadening profile of the laser transition 

POO 

Z(ei,2) = 2i e-"'-2''"-''«i-2"rfx , (18) 



where ^1^2 = (1^1,2 — (^o)/^d is the detuning, normalized to the Doppler width, from the 
transition center for the beams 1 and 2. The independent variables ^1^2 are in correlation with 
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temperature and pressure inside the cavity. For an active medium composed by a gas mixture 
of two isotopes, one has to account for C, and rj and rj', To and F'^,. Here the unprimed 
and primed symbols refer to the A^e isotopes 20 and 22, respectively. In the Doppler limit 
?7 ^ 1 and rj' <^ 1, which is common for middle or large size He-Ne rings, Z{C,) is usually 
approximated as 

ZniO ^ -2ee-«' , 

where the pedices I and R stands for imaginary and real part, respectively. Within the 
above approximations the Lamb coefficients have the following expressions 



ai,2 

0"1.2 



G 



ZiiO) 

=ai,2 + 1^1,2 
fo G 



e 



2 Zj{0) 
TG 



12 



ri2 



Zj{0) 

r/o 



k 



G 



2 Zj{Q) 



[kZn{^i,2) + 

Zl{^l,2) 

Zl{^l,2)U/V 

1 + iUv)' 



k'Zni^l,) 

Zii^'1,2) 



(19) 



k'- 



M^i,2)Uri 

^ 1 + {iLhY . 



where G is the laser gain, F = (7a + 76)/(27a6), = (6 + 6)/2, C = (^1 +^2)/2, Ati,2 are 
the mirror losses experienced by each beam, k and k' are the fractional amount of isotopes 
in the gas mixture. The equation for ^21 and T2\ are obtained from the expression of ^12 and 
ri2 by permuting the subscripts 1 and 2. Substituting ^12 — 0.47, 2 — ~0-49, k = k' = 0.5, 
G ~ 3 X 10~^, rj ~ 0.25, and ^1^2 — 1,48 x 10~^ in Eqs. f|T9|) we get for the closed loop 
operation of G-Pisa (3i — (32 — 10"^"^, 6'i2,2i — 8 x 10~^, ri2,2i — 10~\ and the values in Table 
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